An analytical study of resonant transport of Bose— Einstein condensates 
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We study the stationary nonlinear Schrodinger equation, or Gross-Pitaevskii equation, for a one- 
dimensional finite square well potential. By neglecting the mean-field interaction outside the po- 
tential well it is possible to discuss the transport properties of the system analytically in terms of 
ingoing and outgoing waves. Resonances and bound states are obtained analytically. The transmit- 
ted flux shows a bistable behaviour. Novel crossing scenarios of eigenstates similar to beak-to-beak 
structures are observed for a repulsive mean-field interaction. It is proven that resonances transform 
to bound states due to an attractive nonlinearity and vice versa for a repulsive nonlinearity, and 
the critical nonlinearity for the transformation is calculated analytically. The bound state wave- 
functions of the system satisfy an oscillation theorem as in the case of linear quantum mechanics. 
Furthermore, the implications of the eigenstates on the dymamics of the system are discussed. 

PACS numbers: 03.75.Dg, 03.75.Kk, 42.65.Pc 



I. INTRODUCTION 

Due to the experimental progress in the field of Bose- 
Einstein condensates (BECs), e.g. atom-chip experi- 
ments, achieved in the past few years 0,0,0- it is now 
possible to study the influence of an interatomic inter- 
action on transport. At low temperatures, BECs can be 
described by the nonlinear Schrodinger equation (NLSE) 
or Gross- Pitaevskii Equation (GPE) 



dt 



2m 



V(r)+g\i;(r,t)\ 2 )^v,t) 



(1) 

in a mean- field approach 0, la, S ■ Another important 
application of the NLSE is the propagation of electro- 
magnetic waves in nonlinear media (see, e.g., ||, ch. 8). 
In the case of vanishing interaction (g = 0), Eq. (Q re- 
duces to the linear Schrodinger equation of single particle 
quantum mechanics. 

The treatment of transport within this mean-field the- 
ory reveals new interesting phenomena which arise from 
the nonlinearity of the equation. For example a bistable 
behaviour of the flux transmitted through a double bar- 
rier potential has been found 0. For the same system 
transitions from resonances to bound states can occur 
due to the mean-field interaction 0,0]. However, pre- 
vious studies of nonlinear transport and nonlinear reso- 
nances usually rely strongly on numerical solutions, such 
as complex scaling. There are only few model systems 
of the NLSE which have been studied analytically and 
most of these treatments concentrate on bound states. 
Among these few examples are the infinite square well 
0, 0, 13, the finite square well 0, 0, the poten- 
f, delta-potentials 0,0] and the delta-comb 
an example of a periodic potential. Therefore 
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analytical solutions of the NLSE for other simple model 
potentials arc of fundamental interest. 



In the present paper we solve the time-independent 
NLSE for a simple one-dimensional square-well poten- 
tial in order to analytically indicate and discuss the 
nonlinear phenomena described above. To unambigu- 
ously define ingoing and outgoing waves and thus a 
transmission coefficient |T| 2 , we neglect the mean-ficld- 
interaction g\ip(x)\ 2 outside the potential well. Since 
the one-dimensional mean-field coupling constant g is in- 
versely proportional to the square of the system's radial 
extension aj_ this can be achieved in an atom-chip ex- 
periment by a weaker radial confinement of the conden- 
sate outside the potential well 0, 0, . Especially for 
bound states or very stable resonances our approxima- 
tion is well justified since in these cases the condensate 
density | -0 (a^) | 2 is much higher inside the potential well 
than outside. 



The paper is organized as follows: In section ^ we 
briefly review some basic results for the scattering states 
of the linear Schrodinger equation with a square well po- 
tential which are treated in various textbooks on quan- 
tum mechanics 0, |^, . In section lllll we will then 
analyze stationary scattering states of the NLSE with a 
finite square- well-potential. In section lTvl we again have a 
look at the linear system with respect to bound states in 
order to prepare for discussing the transition from reso- 
nances to bound states for an attractive (section[VJ and a 
repulsive mean-field interaction (section lVIfl . In addition 
to our discussion in terms of stationary states we have a 
brief look at the dynamics of the system in section IV11I 
where we solve the time-dependent NLSE numerically 
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II. LINEAR SCHRODINGER EQUATION: 
SCATTERING 



The scattering states of the linear time-independent 
Schrodinger Equation 



ti 2 d 2 

2m dx 2 ^ 



tp(x) = Eip(x) 



(2) 



for a particle with mass m where the potential is given 
by 



V(x) 



V , \x\ < a 
0, \x\ > a 



(3) 



with Vb < arc obtained for positive particle energies 
£7 > 0. We assume an incoming flux from x = — oo and 
make the ansatz 



Ae ikx +Bc -ikx x < _ a 

Fe iqx + Ge~ iqx \x\ < a 
C e lkx x > a , 



tp(x) 



(4) 



where the wavenumbers inside (|a;| < a) and outside 
(| a; | > a) the potential well are given by 



- V )/n. 



and 



k = V2mE/h, 



(5) 



(6) 



respectively. The amplitudes of transmission and reflex- 
ion are then determined by 



T := C/A and R := B/A 



(7) 



They satisfy the relation |i?| 2 + |T| 2 = 1 as the 
Schrodinger equation is flux-preserving. We have to make 
the wavefunction and its derivative continuous at x = ±a. 
Together with the normalisation condition A = 1 these 
boundary conditions determine the wavefunction of the 
system. By eliminating F and G we obtain the amplitude 
of transmission 



T(E) = 



1 



exp(— 2ifca) 



cos(2ga) 1 - (i/2)[(g/fc) + (k/q)] tan 2 (ga) 
and thus the transmission probability 

|W=(l + i(|-J) 2 sin^a)) 1 



(8) 



(9) 



We define a resonance as a maximum of the transmission 
coefficient where the potential well is fully transparent, 
that is \T(E)\ 2 = 1. Then Eq. © implies the resonance 
condition 
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FIG. 1: (Color online) Transmission probability (left) and 
scattering phase (right) for the potential Vb = — 10 and a = 2 
(for h = m — 1). 



The resonance energies are then given by 



H 2 q 2 H 2 ir 2 , 



2m 



8ma 



(11) 



where n is an integer number and sufficiently large to 
make E n > 0. A Taylor expansion of Eq. © in terms of 
E — E n yields 



\T{Et 



where T n is given by 



(r„/2) 



(E - E n ) 2 + (r„/2)2 



2\JE n (E n — Vb) 
2E n - V ' 



(12) 



(13) 



This means that the transmission coefficient can be ap- 
proximated by a Lorentzian with width T n in the vicinity 
of the resonances. We define the scattering phase 9(E) 
by 



2qa = WK. 



(10) 



T{E) = \T(E)\c 



W(E)-2ika 



(14) 
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and obtain 



tan 0(E) = -(- + -) tan(2ga). 
2 \k qj 

In the vicinity of a resonance this reduces to 



6(E) w arctan ( — (£ - E n ) 

-L n 



(15) 



(16) 



Note that the scattering phase is zero for E = E„ . The 
bound states of this system will be discussed in section 

IVl 



III. NLSE: SCATTERING 

Now we consider the same square well potential with an 
additional mean-field term <7|^(x)| 2 inside the potential 
well so that for |x| > a the wavefunction must satisfy the 
free linear Schrodinger equation 



2m dx 2 



(17) 



and for |x| < a the time-independent NLSE 
2m dx 2 



g\^{x)\ J +V -fi)i>{x) = (18) 



with a constant potential Vq < 0, where fi is the chemical 
potential. 

For an incoming flux from x = — oo we now make the 
ansatz 



ip(x) = 




(19) 



with the wavenumber JfjJ). Inside the potential we assume 
an intrinsically complex solution 



f(x) = y/S&j*™ 



(20) 



where S(x) and are real functions. As already 

shown by Carr, Clark and Reinhard [T3. ll3T | a complex 
solution of the NLSE with a constant potential is given 
by the equations 



&(x) 



a 



and 



S(x) = e + ip dn 2 (gx + 5\p) 



(21) 



(22) 



in terms of the Jacobi elliptic function dn with the real 
period 2K(p), where K(p) is the complete elliptic int egra l 
of the first kind and p is the elliptic parameter [2fill27lf28l |. 

The real parameters a,p, e, ip, g and S have to satisfy the 
relations 



gm 



(23) 



3 h 2 . . 2 



(24) 



-— (q 2 V e( P - 1) - a 2 ) + ge 3 -fa- V )e 2 = . (25) 
2m v ' 

Again, the wavefunction and its derivative must be con- 
tinuous at x = ±a. At x = +a we obtain the equations 



/(o) = Cc ite , f'(a) =ikCe ika 
By inserting Eq. (121)1) we arrive at 



— S (a) + ia = ifc S(a) . 



(26) 



(27) 



Since we are considering positive chemical potentials 
(i > 0, the wavenumber k is real so that we obtain the 
conditions 



and 

At x — — a we obtain 

^—iki 



S'(a) = 



a = kS(a). 



and 



ik(Ai 



-ika 



Be 



ika\ 



/'(-«)• 



(28) 



(29) 



(30) 



(31) 



With the equations (|2U|I and l|29|) these two conditions 
can be written as 



y/S(-a)Ae 



-ika 



i 



S(-a) + S(a) - i 



S'(-a) 
2k 



(32) 



y/S(-a)Bc 



-\-ika 



S(-a) - S(a) +i 



.S'(-a) 



2k 



3 i*(-o)_ 

(33) 

As in the linear case, we define the amplitudes of re- 
flection and transmission as R := B/A and T :— C/A. 
Equations (j32() and (|33|) then lead to the reflection prob- 
ability 

, R|2 = [S(-a)-S(a)} 2 + S'(-a) 2 /(4k 2 ) 

11 [S{-a) + S(a)] 2 + S'{-a) 2 /(4k 2 ) [ ' 

and the transmission probability 

45(o)5(-a) S(a 



\T\ 



[S(-a) + S(a)} 2 + S"(-a) 2 /(4fc 2 



\A\< 



(35) 
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It can easily be verified that \R\ 2 + \T\ 2 = 1. The scat- 
tering phase 9(n), defined by T(p) = |T(^i)|cxp(i0(/z) — 
2ifca), can be expressed as 



0(H) = arg(C) - arg(A) + 2kc 



(36) 



+Q kS(a) , S'(-a) 
ax + arctan - 



-a S(x) 



2k(S(-a) + S(a)) 



The condition 1)28|) can be further exploited by rewrit- 
ing it as 

—2ipgp dn(u\p) sn(u\p) cn(u|p) = , (37) 

with the abbreviation u = ga + S. As the dn-function is 
always nonzero, the solutions of Eq. are given by the 
zeros of either of the two functions cn(u\p) or sn(u\p). 



Case 1 : sn(ga + 8\p) = 

^ga + S = (2j + l)K(p) 

=>S={2j + l)K(p)-ga. 



(38) 



(39) 



The variable j denotes an integer number. This 
means 

dn 2 (ga + 5\p) = dn 2 (2jK(p)\p) = 1 (40) 



or 



S(a) =e + (p. 

Case 2 : cn(ga + S\p) = 

In analogy to case 1 we get 



and 



S = (2j + l)K(p) - ga 



S(a)=£ + cp(l-p). 



(41) 



(42) 



(43) 



For given values of the chemical potential /z and the 
incoming amplitude \A\, the parameters p,e, and tp = 
— h 2 g 2 /gm are completely determined by the equations 
iJUJl and {231 an( l the absolute square of Eq. {32$ . 

16k 2 \A\ 2 S(-a) = Ak 2 [S(-a) + S(a)} 2 + S'(-a) 2 . (44) 

The positions of the resonances arc defined by \T\ 2 = 1 
or |i?| 2 = 0. With Eq. this is equivalent to the 
conditions 



and 



S'(-a) = 



S(-a) = S(a) . 



(45) 



(46) 



Now we consider these conditions for the two cases 
mentioned above. 



Case 1 : sn(ga + 6\p) = or S = (2j + l)K(p) — ga 
The condition (|45f) now reads 

S (—a) = —2ipgpcn(u\p) sn(u\p) dn(u\p) = (47) 

with u = —ga + 6. The solutions of Eq. (|45|l arc 
again given by the zeros of either the sn-function 
or the cn-function. One can show easily that only 
the zeros of the sn-function are compatible with the 
condition Thus we obtain 



-2ga + 2jK(p) = £K(p) 



(48) 



where I and j are integer numbers. The resonance 
positions arc then given by 



K(p) 
g = n, 



(49) 



where the integer number n is defined as n = £ — j . 

Case 2 : cn(ga + S\p) = or 6 — 2jK(p) — ga 

In analogy to case 1 one can show that the reso- 
nances are also determined by Eq. I)49|l . 

In the case of a resonance there is no reflected wave 
(B = 0) so that the equations (SJ, @3), I® and (IBjl 
lead to 



\A\ 2 = S(-a) = S(a) = 



e + ip 



case 1 



e + ip(l — p) case 2 



or 



e(p) 



\A\ 2 



case 1 



\A\ 2 - ip(l - p) case 2 



(50) 



(51) 



With the help of the equations lj5T)l . and $Zty we ob- 
tain a formula for the chemical potential of a resonance, 



/ n „ 3 , ,,2 ft2 k2 (p) 2 
M P ) = V Q + - 9 \A\ 2 + --^>n 2 . 



(1 + p) case 1 
(1 — 2p) case 2 . 

(52) 

The chemical potential depends on the parameter p and 
the integer number n which must be sufficiently large to 
make fiR ^ 0- As in the linear case (cf. Eq. I|ll[l ) there 
is a term which grows quadratically in n. This term can 
be viewed as an effective increase of the chemical poten- 
tial in comparison to the energy of the linear system in 
case 1 and respectively as a decrease in case 2. How- 
ever, it turns out (cf. below) that the energy shift due to 
nonlinear interaction is mainly determined by the term 
|g|^4| 2 which is a direct consequence of the effective in- 
crease or decrease of the potential well by the mean-field 
term <7|^(a;) | 2 . As we will show below, we have p — > if 
g — > so that K(p) — > 7r/2 and Eq. (|52[) reduces to the 
equivalent equation 111|) for the resonance energies of the 
linear system. 

The elliptic parameter p can be determined by insert- 
ing the expressions l|52|l and 151|l into Eq. (|25|l . We ob- 
tain: 
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FIG. 2: (Color online) Squared modulus of the wavefunctions 
of the two most stable resonances of the potential Vo = — 10, 
a — 2 for the attractive nonlinearity g = — 2. The normaliza- 
tion \A\ = 1 has been used. 



FIG. 3: (Color online) Squared modulus of the wavefunctions 
of the two most stable resonances of the potential Vo = —12, 
a = 4 for the repulsive nonlinearity g — +5. The normaliza- 
tion \A\ = 1 has been used. 



Case 1: 



pK 4 (p)=g\A\ 2 (V +g\A\ 2 ) 



2m 2 a 4 
h 4 n 4 



(53) 



For g = this equation can only be solved by p = 0. 
The inequality pK 4 (p) > implies 



g\A\ 2 (V + g\A\ 2 )>0. 



(54) 



Since Vo is negative this means either g < or 
g\A\ 2 > | Vb|. For small values of p the left hand 
side of (|53|) can be approximated by 



P K\p)=(^f (p + p 2 ) + 0(p*) 



1 2g\A\ 2 (V +g\A\ 2 )m 2 a* 
4 (ir/2) 4 h 4 n 4 




(55) 



(56) 



Case 2: 



Otto 

p{l-p)K 4 (p)^- g \A\ 2 (V +g\A\ 2 ) 1 ^- r . (57) 



For g = 0, this equation is solved by p = 0. The 
solution p = 1 is incompatible with Eq. be- 
cause it leads to a negative chemical potential [ar. 
From p(l - p)K 4 (p) < we obtain the condition 
< g\A\ 2 < \Vq\. For small values of p, the left 
hand side of l|57(l can be approximated by 



p(l-p)K 4 {p) 



p + 0(p 3 ) 



and we get 



2 g \A\ 2 (Vo+g\A\ 2 )m 2 a 4 
(Tr/2) 4 h 4 n 4 



(58) 



(59) 



Case 2 is valid for < g < |Vb|/|A| 2 and case 1 is valid 
for any other value of g. Obviously the parameter p and 
the chemical potential pr only depend on the product 
rjlAj 2 . Without loss of generality we henceforth assume 
\A\ 2 = 1 in all figures and numerical calculations. The 
squared modulus of the resonance wavefunctions inside 
the potential well \x\ < a in dependence of p and n is 
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given by 
\M')\ 2 = \A\ a -- 8 



[dn 2 (gx - nK(p) \p) - l] (60) 



in case 1 and 



h 2 g 2 



\^ R (x)\ 2 = \A\ 2 [dn 2 {gx + (1 - n)K(p)\p) 

9 



- 1+p] 

in case 2 with g = nK(p)/a. 



(61) 




FIG. 4: (Color online) Transmission coefficients for different 
attractive nonlinearities for the square well potential Vq = 
— 10 , a = 2. The positions of the resonances, calculated with 
the formula 15211 . are marked by '*'. 

Figures |21 and [3] show the squared modulus of reso- 
nance wavefunctions for attractive respectively repulsive 
nonlinearities in the region of the potential well \x\ < a 
where scaled units with h = m = 1 are used, as for all 
figures and numerical calculations in this paper. The 
functions are symmetric with respect to x = and have 
no nodes. The quantum number n indicates the number 
of the minima. The smaller the chemical potential ^.r(ti) 
is, the smaller is the value of |^(x)| 2 & t its minima. 

Figures 01 and [3] show the transmission coefficient 
|T(/i)| 2 in dependence of the chemical potential /z for the 




FIG. 5: (Color online) Transmission coefficients for different 
attractive nonlinearities for the square well potential Vo = 
— 10 , a = 2. The positions of the resonances, calculated with 
the formula 152H . are marked by '*'. 



potential Vq = —10 and a — 2 for different attractive 
nonlinearities g < 0. For small values of \g\, the be- 
haviour of the transmission coefficient is very similar to 
the linear system (see Fig. ^1 . With increasing attractive 
interaction, the curves bend more and more to the left 
so that the resonances are shifted towards smaller val- 
ues of (i. This shift is mainly determined by the term 
|g|^4| 2 in Eq. For g = —2, the lowest resonance so 

far has disappeared into the bound state regime fx < 0. 
Between g = — 6 and g = — 8 the next resonance dis- 
appears. It can be shown that these resonances become 
bound states of the system (see section Wl . For higher 
values of \g\ the bending of the curve in the vicinity of 
the resonances leads to a bistable behaviour. For a fixed 
value of /x, there are now different wavefunctions which 
lead to different transmission coefficients \T(fi)\ 2 . Due 
to the nonlinear mean field interaction the state adopted 
by the system for an incoming flux from x = — oo is not 
only determined by the chemical potential \i and the am- 
plitude A but also by the density distribution \i/j(x,t)\ 2 
of the condensate inside the potential well so that the 
system has a memory (hysteresis). This will be further 





°0 2 4 6 8 10 

FIG. 6: (Color online) Transmission coefficients for different 
repulsive nonlinearities for the square well potential Vo = — 12 
, a = 4. The positions of the resonances, calculated with the 
formula 152H . are marked by '*'. 



2 4 6 8 10 

H 

FIG. 7: (Color online) Transmission coefficients for different 
repulsive nonlinearities for the square well potential Vo = —12 
, a — 4. The positions of the resonances, calculated with the 
formula 152H . are marked by '*'. 



investigated in section IVlTl 

Now we consider repulsive nonlinearities g > 0. Fig- 
ures [5] and [7] show the transmission coefficient \T\ 2 in 
dependence of /j, for a potential Vb = — 12, a = 4 for 
different repulsive nonlinearities. Here we observe a shift 
of the resonance positions towards higher values of fi, 
since the term §<?|^4| 2 is now positive. New resonances 
appear for increasing \g\ which originate from the bound 
states of the linear system (see section HVj) . Now the 
curves bend to the right and, as in the attractive case, 
the transmission coefficient shows a bistable behaviour. 
Between g = 4.5 and g = 5, we observe a bifurcation 
phenomenon at fi as 1 which resembles an avoided cross- 
ing. This turns out to be a beak-to-bcak scenario which 
also occurs in two-level systems. These phenomena will 
be discussed in detail in a subsequent paper |29l |. In all 
transmission curves we observe that, for a fixed value of 
<7|A| 2 , the effect of the nonlinearity on the shape of the 
resonances is stronger for small values of /x than for large 
ones, because the kinetic energy is higher and the mean- 
field energy has a comparatively smaller influence. This 
also manifests itself in formula for [1r where the 



term proportional to n 2 counteracts the dominant term 

Figures ISlllI show scattering phases in dependence of 
the chemical potential fj, which have been calculated by 
numerical integration of Eq. (|37|l for the parameters of 
the figures 14171 For weak nonlinearities the scattering 
phase behaves similar to the scattering phase of the linear 
system which is zero at the positions fi = fiji of the res- 
onances. For stronger nonlinearities the resonance scat- 
tering phase increasingly deviates from zero. Since the 
squared modulus of the wavefunction enters Eq. I|37[) . the 
scattering phase 9(fi) inherits the bistable behaviour of 
the transmission coefficient |T(^)| 2 . Wherever beak-to- 
beak structures occur in the transmission coefficient the 
scattering phase shows loops. 

IV. LINEAR CASE: BOUND STATES 

Now we consider bound states of the square well poten- 
tial. Again we first have a brief look at the linear system. 
Bound states are found for energies Vq < E < 0. The 
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FIG. 8: (Color online) Scattering phases for different attrac- 
tive nonlinearities for the potential Vo — — 10, a = 2. The 
positions of the resonances, calculated with the formula 1521 . 
are marked by dotted vertical lines. 



area |ar| > a outside the potential well is classically for- 
bidden and the wavefunction decays exponentially with 
the coefficient 



K = V-2mE/h. (62) 

Inside the well it still oscillates with the wavenumber 

q = ^2m(E- V )/h. (63) 

Due to the symmetry V(x) = V(—x) of the potential we 
can choose eigenfunctions which are either of odd or even 
parity. For even parity we make the ansatz 



tp + (x) ■■ 
and for odd parity 

i>-(x) = 



7+ e KX x < —a 
■q + cos(qx) \x\ < a 
7+ c~ KX x > a 



7_ e KX x < —a 
T]- sin(ga;) \x\ < a 

— 7_ e~ KX x > a 



(64) 



(65) 
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FIG. 9: (Color online) Scattering phases for different attrac- 
tive nonlinearities for the potential Vo = — 10, a = 2. The 
positions of the resonances, calculated with the formula l|52lL 
are marked by dotted vertical lines. 



As in the case of scattering states the wavefunction and 
its derivative must be continuous at x = ±a. For even 
parity we get the condition 



qtan(qa) = k, 



(66) 



which determines the energy of the bound states. For 
odd parity we obtain 



qcot(qa) = — k . 



(67) 



The transcendental equations 16611 a nd (|67|l can be 
solved numerically. One can show [24| that the num- 
ber of symmetric or antisymmetric solutions is given by 



and N- 



Rq 
7T 



(68) 



respectively, where [y] > denotes the smallest integer 
greater than y and i?o is given by 



(69) 
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FIG. 10: (Color online) Scattering phases for different repul- 
sive nonlinearities for the potential Vo = — 12, a = 4. The 
positions of the resonances, calculated with the formula Q52I . 
are marked by dotted vertical lines. 



FIG. 11: (Color online) Scattering phases for different repul- 
sive nonlinearities for the potential Vo = — 12, a = 4. The 
positions of the resonances, calculated with the formula l|52lL 
are marked by dotted vertical lines. 



V. ATTRACTIVE NONLINEARITY: BOUND 
STATES 

In this section, we consider bound states of the non- 
linear system described by the equations (|18|) and 117|) 
for attractive nonlinearities g < 0. In analogy to the lin- 
ear case, bound states occur only for chemical potentials 
H < 0. As in the linear case, we look for solutions of even 
and odd parity. We make the ansatz 



7 + e KX x < —a 
">P( X )+ = { i +™(q + x\p + ) \x\ < a 



7+e" 



x > a 



for even and 



7_ e KX x < —a 

ip(x)- = ( I-Cii(q-X + K (p_)|p_) \x\ <a 



-7_e 



x > a 



for odd parity with 



(70) 



(71) 



The solutions inside the potential well \x\ < a are given 
by Jacobi elliptic functions cn with amplitudes 



/+ = W — 



'Q±P± 
gm 



(73) 



and wavenumbcrs 



Q± 



I 2m(fi - V) 

h 2 (i-2 P± )- 



(74) 



1. symmetric solutions 

The continuity of the wavefunction and its deriva- 
tive at x = ±a yields with the abbreviation u = 
Q+a) 



7+e Ka = 1+ cn(w|p + ) 



(75) 



k = \J — 2mfi/h. 



(72) 



«7+e 



I + q + sn(u|p+)dn(u|p+). (76) 
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By inserting Eq. (|75|) into Eq. l(7fi|) we get the con- 
dition 

Kcn(u|p+) — q+ sn(u|p + )dn(u|p + ) = (77) 

for symmetric bound states. In the limit p — > 0, the 
Jacobi elliptic functions converge to trigonometric 
functions and we regain the condition (|66(l of the 
linear system. 

2. antisymmetric solutions 

In analogy to the case of even parity we obtain the 
condition 

«;cn(u|p_) — Q_ sn(u|p_)dn(it|p_) = (78) 

for antisymmetric bound states with u = g_ a + 
K(p_), which also converges to Eq. (|67(l in the limit 
p^ 0. 

Now we show explicitly that the resonances discussed 
in section HTTI become bound states for a sufficiently strong 
attractive mean-field interaction. Due to the behaviour 
of the wavefunction outside the potential well, transitions 
from resonances to bound states occur at (j, = 0. For 
H = the condition (|77(l for symmetric bound states 
becomes 

sn(g + a\p+) = (79) 

which is solved by 

g+ =2j+ K(p+)/a (80) 
with integer j+ . Inserting Eq. (j80(l into l|74|) yields 

m\V \a 2 



{l-2p + )K\p + ) = 



2h 2 j% ■ 



(81) 



For small values of p+ we can make a Taylor approxima- 
tion and obtain 




§) +C + + o(pD 



(82) 



as an approximate solution for the parameter p + where 
the transition occurs. Here we have introduced the ab- 
breviation 



•32 
23 



1 



2m\V \a 2 

H 2 TT 2 jl 



32 
23 



1 - 



3% 



(83) 



The quantity Rq /tt in £ + determines the number of states 
N + = [Rq /tt] > for .g = (see sction IIV|I . Since p + must 
be positive, C+ must be positive as well and we obtain 
the condition 



3+ > N + 



(84) 



Next we want to demonstrate that the wavefunction of 
a symmetric bound state continuously merges into a res- 
onance wavefunction with an even quantum number n at 



fi = if we use an appropriate normalization. To avoid 
confusion, all the parameters of the resonance wavefunc- 
tion now carry an additional index R. For fi = we 
have k = k = (see equations © and so that the 

squares of the wavefunctions outside the potential well 
are identical if we use the normalization 



17+1 = 

The equations l|75)l. (JTSJ) and (jSOJl then imply 
T |2 Ah 2 j 2 + K(p + ) 2 P+ 



\A\< 



gma A 



(85) 



(86) 



Inside the potential well the squared modulus of the res- 
onance wavefunction is given by 



l^(z)| 2 = \Af 



h 2 K 2 {p R )n 2 



gma^ 



nK{pi 



dn I x\p R J - 1 

(87)" 

(see Eq. (|60|) ) as g < corresponds to case 1 in section 
IIIII Because of k = 0, the phase of the resonance wave- 
function is constant and the phases of ip R and ip+ can be 
chosen identical. Using Eq. Q86|l and the relation between 
the squares of the Jacobi elliptic functions [26l I27I |28| we 
arrive at 



l<M*)l 2 = - 



hHj 2 + K(p + ) 2 P+ h 2 K 2 (p R )n 2 



gma^ 

2 / nK(p R ) | 
cn x\p R 



gma^ 
- 1 . 



-PR 



(88) 



By comparing \i/j r (x)\ 2 with the squared modulus of the 
bound state wavefunction, 



liMs)| s 



h 2 AjlK(p^ 



'P+ 2 
cn 



gma* \ a 

we see that \ijj R {x)\ 2 is equal to |?/> + (:z;)| 2 if 
PR =P+ ='-PT+ 

and 

n = 2j+ , 



-x\p+ 



(89) 
(90) 

(91) 



where the index T denotes transition. Eq. i|86|l implies 
that the transition occurs for the nonlinear interaction 
strength 



9T+ 



h 2 n 2 K{p T +) 2 PT+ 
ma 2 \A\ 2 



(92) 



Inserting Eq. {H21> and (JHTJ into Eq. for [i R and 

Eq. H53JI ■ which connects g and p R , we see that 



Vr(P = PT+) = 



(93) 



and that Eq. Ij53(l holds. Thus we have proved that 
Eq. I|90|) is valid so that resonances of even quantum num- 
ber n indeed become symmetric bound states. 
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symmetric (n 
numerically exact 


= 6) 
approx. 


antisymmetric ( 
numerically exact 


re = 7) 
approx. 


PT 


0.0644 


0.0643 


0.2044 


0.205 


9T 


-1.4772 


-1.4749 


-6.9159 


-6.969 


9t 


-0.7325 


-0.7331 


-3.3592 


-3.346 



TABLE I: Critical parameters for the transition of the two 
most stable resonances of the potential Vo = — 10, a — 2 to 
bound states. The approximations 182H for p+ and 19711 for 
p- have been used. 




0.8 



0.6 



0.4 



0.2 



n=7 



-2 



FIG. 12: (Color online) Squared modulus of the wave- 
functions of the two most stable resonances of the potential 
Vo = —10, a = 2 in each case exactly at the critical points for 
the transition from resonances to bound states (comp. table 
HJ. Here, the minima are also zeros of the wavefunction (cf. 
Fig. |2J . The normalization \A\ = 1 has been used. 

In a similar way one can show that resonances of odd 
quantum number n become antisymmetric bound states 
for a sufficiently strong attractive nonlincarity g < 0. We 
obtain (cf. Eq. JEU) 

Q- = (2j- + l)K(p-)/a (94) 
with integer j_ and n = 2j_ + 1 . If we choose the 



normalization |7_| = \A\, the transition occurs at 

h 2 n 2 K{p T _) 2 p T _ 
^ T ma 2 \A\ 2 

The parameter pt- '■= P- = Pr is determined by 



(95) 



f-^^'-s^? (96) 



Pt- 



where we have introduced the abbreviation 



C- 



32 
23 



8m\V \a 2 



ft 2 7r 2 (2j_ + iy 



(98) 



The condition p?- > implies £_ > and thus j_ > iV_ 
where = [-^f 1 — |] is the number of antisymmetric 
bound states for g = as defined in Eq. (|68|l . 



0rr= 



-0.1 - 



-0.2 



-0.3 



-0.4 



-0.5 



-0.6 



-0.7L 



0.1 



0.2 
P 



0.3 



0.4 



FIG. 13: Dependence of the effective nonlinear interaction 
strength on the parameter p± for the symmetric bound 
state n = (solid line) and the antisymmetric bound state 
71 = 1 (dashed line). 

Figure lT2l shows the squared modulus of wavefunctions 
at the transition from resonances to bound states. As in 
the case of resonances the functions arc symmetric and 
have n minima (cf. Fig. |5J). The quantum number n 
counts all the resonances and bound states defined above 
and is equal to the number of minima of the squared 
modulus of the wavefunction. For bound states, these 
minima are also nodes of the wavefunction. Thus, the 
famous oscillation theorem of linear quantum mechanics 
(cf. jlSl) is a l so valid in the nonlinear case considered 
here. 

In order to study the shift of the bound state energies 
due to mean-field interaction it is useful to introduce an 
effective nonlinear interaction strength 



(99) 
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-2 
0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 











n=1 






g eff =-0.15 











FIG. 14: (Color online) Wavefunctions of the two lowest 
bound states of the potential Vo = —10, a = 2 for g cB = —0.1 
('-') in comparison with the respective wavefunctions of the 
linear system ('- -'). The normalization J® a \ip(x)\ 2 dx = 1 has 
been used. 



(cf. [lSj ) in the interaction region \x\ < a. The integra- 
tion can be carried out analytically. We obtain 



-off 



9- 



^[E{ e -a + K(p-)\p-)-E(p-) 
ma 

(i- P _)e-o] (loo) 



for antisymmetric and 
h 2 g+ 



9f 



[E(g + a\p + )-(l-p+)g + a] (101) 



for symmetric bound states where E(p) and E(u\p) are 
the complete and incomplete elliptic integral of the sec- 
ond kind, respectively. The critical effective nonlinear 
interaction strength for the transition from resonances to 
bound states is given by 



off 

9t± 



9t ± \A? 

PT ± 



E( PT± ) 
K (PT ± ) 



(1-PT ± ) 



(102) 



As seen from Eqs. <|100[) and (|101|) . \g (p±)| increases 
faster for bigger wavenumbers g± (respectively bigger 



=t -6. 




-12 



FIG. 15: (Color online) Shift of the bound states towards 
lower chemical potentials in dependence of |g eff | for p± < 0.99 
for the potential Vo = — 10, a = 2, for symmetric (x), and 
antisymmetric (o) bound states (attractive nonlinearity) . 



chemical potentials n, cf. Fig. I13|) . This means that for 
a given value of g , states which are strongly bound are 
more affected by the nonlinear interaction than weakly 
bound states. As in the case of resonances this is due 
to the fact that for weakly bound states the energy of 
the mean-field interaction is much smaller than the ki- 
netic energy and thus has a comparatively weak influ- 
ence. This can also be observed in Fig. ^] which com- 
pares the wavefunctions of the two lowest bound stats of 
the potential V = -10, a = 2 for g cS = -0.15 with the 
respective wavefunctions of the linear system (g eff = 0). 
The wavefunction of the state n = differs much more 
clearly from its linear counterpart than the wavefunction 
of the state n = 1. 

Figure lT31 shows the shift of bound states in dependence 
of \g \- The appearance of new bound states can be 
observed. For large values of |g cff |, which correspond to 
elliptic parameters p > 0.5, the states all lie below the 
potential well depth Vo = —10- For small values of p± the 
chemical potential fj,(g±) can be approximated linearly 
in terms of g*± , 



h 2 g 2 _p. 
2ma 



cn 2 (g_x + K(p_))dx (103) 



l'l-^)+0(p?.), 



where 



Eo = V + h 2 q 2 /(2m) 



(104) 



is the energy of the respective bound state of the linear 
system. For symmetric states we get 



gf = -p + (Eo-V )[l 



sin 



(2ga) \ 



2qa J 



0<j>l). (105) 
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FIG. 16: (Color online) Shift of the bound states towards 
lower chemical potentials in dependence of \g eB \ for p± < 0.4 
for the potential Vo = —10, a — 2. 'o': numerically exact 
values respective linear approximation due to Eq. 11091 
(attractive nonlinearity) . The sudden disappearance of the 
curves at g cB (p = 0.4) is due to the fact that in this figure the 
Jacobi elliptic parameter was "artificially" limited to p < 0.4 
in order to avoid confusion and for a better comparison of the 
exact values and the approximation. 



Then, for small values of p± , we have 

sin(2ga) \ 



P± 



-9f 



(E - V ) 1 ± 



2qa J 

and the chemical potential is given by 



2m 



(l-2p±), 



(106) 



(107) 



where g± is the wavenumber of a Jacobi elliptic function 
en with quarter period K(p±). In the limit p± —> 
the cn-function becomes a trigonometric function with 
quarter period ir/2 and wavenumber q. This motivates 
the approximation g± w 2qK(p±)/% which leads to 



h 2 q 2 4 
2m 7T 2 



(l-p±)K 2 (p±). 



(108) 



Inserting the equations (1 104(1 . ((106j) and expanding 
K 2 {p)(l - 2p) = (tt/2) 2 (1 - |p) + C(p 2 ) we obtain 



1± 



sin(2ga) 
2qa 



(109) 



A similar relation was derived for a delta-shell potential 
[l8j . For small values of p± this approximation is in good 
agreement with the numerically exact calculations (cf. fig 
I16fl . In the limit Vo — » — oo and thus q — ► oo we obtain 
the formula 



(110) 



for an infinitely deep square well potential 0, fl3L Il3| . 



VI. REPULSIVE NONLINEARITY: BOUND 
STATES 

As in the previous section we look for solutions of odd 
or even parity. We make the ansatz 



7-1- e KX x < —a 

"4>(x)+ = { 1+ sn(p + x + K{p + )\p+) \x\ < a 



(111) 



7+e 



x > a 



for even solutions and 



il>{x). 



j_e KX x < —a 
Isii(q-x\p-) \x\ < a 
—j_e~ KX x > a 



(112) 



for odd solutions where k = y 7 — 2m/i/H. 

For |x| < a the solutions arc given by Jacobi elliptic 
sn-functions with amplitude 



and wavenumber 



Q± 



"Q±V± 
gm 



I 2m(n-V) 

h 2 {i+p±y 



(113) 



(114) 



Again the wavefunctions and their derivatives must be 
continuous at x = ±a. Thus we obtain the equations 

k sn(it|p+) + g+ CQ.{u\p + )<hi{u\p + ) = (H5) 

for even parity, with u = g+ a + K(p+), and 

«sn(g_ a\p-) + g~ cn(g_ a|p_)dn(p_ <AP-) = (H6) 

for odd parity, which determine the bound states of the 
system. In the limit p — > the conditions (|115fl an d 
(|116fl converge to the equations and l|67|) of the linear 
system. 

For the case of an attractive nonlinearity, we have 
shown that resonances with even/odd quantum number 
n become bound states of even/odd parity due to the at- 
tractive mean-field interaction. In a similar way one can 
show for the case of repulsive nonlinearity that even/odd 
bound states become resonances with even/odd quantum 
number n. 

As before, we use the normalization |7±| = |^4|. Thus 
we obtain the conditions 



where 



g±a = nK(p±) 



2j + , j+ G N even parity 
2j_ + 1, j_ G No odd parity. 



(117) 



(118) 



14 



The transition occurs at 

h 2 n 2 K (pr± ) 2 pr± 
9T± = ml 2 \A\ 2 ' 

where the transition parameters p± are given by 
{l+PT±)K 2 {p T ±) 



h 2 n 2 



(119) 



(120) 




OA 

-) -C(n)+0( PT± ) 3 



with 



32 / 8m\V \a 2 



25 



h 2 TT 2 



(122) 



Since px ± has to be positive this implies the conditions 



j± < N± , 



(123) 



where N + and are the numbers of even and odd 
bound states of the linear system. This expresses the ob- 
vious fact that the number of newly created resonances 
cannot be greater than the number of bound states of the 
linear system. 

As in the case of an attractive nonlinearity, the quan- 
tum number n numbers the the resonances as well as the 
symmetric and antisymmetric bound states and is equal 
to the number of the minima of the squared modulus 
of the wavefunction. In the case of bound states, these 
minima are also zeros of the wavefunction, so that the 
oscillation theorem is valid. 




' eft 



FIG. 17: (Color online) Shift of the bound states towards 
lower chemical potentials in dependence of g eS for p± < 0.995 
for the potential Vb = —10, a — 2. '+': symmetric bound 
states, V: antisymmetric bound states. The sudden disap- 
pearance of the curves at g cfl (p = 0.995) is due to the fact 
that in this figure the Jacobi elliptic parameter had to be 
"artificially" limited to p < 0.995 for numerical reasons. 



In terms of the effective nonlinearity defined by 
Eq. O, we obtain 



9- 



h 2 e . 



[g-a - E{g-a\p-)} 



(124) 



for antisymmetric bound states and 
H 2 g+ 



[q+ a + E{p+) - E{g+ a + K(p + )\p+)] 

(125) 

for symmetric bound states. For the transition we get in 
both cases 



9t ± 



9t ± \A? 

PT ± 



E(PT ± ) 

K(p T± ) 



(126) 



The magnitude of g eS has the same qualitative depen- 
dence on p± as in the attractive case so that the nonlin- 
earity has a stronger effect on lower bound states. 

Figure IT7I shows how the bound states of the system 
are shifted towards higher chemical potentials if g eS is 
increased, until they reach fi = and develop into reso- 
nance states. Approximating the chemical potential lin- 
early in terms of g eS results in 



sin(2ga) \ 1 ff 
2qa J 9 



(127) 



which is identical with the respective Eq. I|109fl for at- 
tractive nonlincaritics. 



VII. DYNAMICS 

In this section we briefly discuss the implications of the 
results for the nonlinear eigenstates presented above on 
the dynamics. To this end we assume that a source out- 
side of the square-well potential emits condensed atoms 
with a fixed chemical potential [i. In fact, we numeri- 
cally solve the time- dependent NLSE with an additional 
delta-localized source term (cf. 0) 



,dtp(x,t) 
1 di 



2 8x 2 

S(t) exp(— i/zt)J(x — xq) , 



+ V(x)+g{x)\^{x,t)\'- 



ip(x,t) 
(128) 



where units with H = m = 1 are used. As above, V(x) = 
V H(a — |x|) is a square- well potential with depth Vo = 
— 12 and half- width a = 4. The nonlinearity is non-zero 
only within the square- well g(x) = gH{a — \x\). Here, 
H{x) denotes the Hcaviside step function. The source 
is located at xo = —15.5. The initial state is an empty 
waveguide tp(x, t = 0) = 0. The source strength S(t) 
is then ramped up adiabatically, so that finally a steady 
transmitting solution tp(x,t en d) is found. To determine 
the incident current and the transmission coefficient, a 
superposition of plane waves with momenta k = ±y/2jl 
is fitted to the wave function ip(x,t en d) for xq < x < 
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—a (upstream region) and x > a (downstream region). 
The NLSE is numerically integrated in real space using a 
predictor-corrector step [23|. To avoid reflections at the 
edges of the computational interval, one-way boundary 
conditions are used as described in |32l |. An example 
of the wave function calculated in this way is shown in 
Fig.[H 




FIG. 18: Transmission of the square-well potential Vo = 
— 12, a = 4 with a repulsive nonlinearity g\A\ 2 = 4. Shown is 
the transmission coefficient of the stationary states calculated 
as described in the previous sections (solid lines) in compari- 
son to the transmission coefficients obtained from a numerical 
solution of the time-dependent NLSE as described in the text 
(crosses) . 

The results of these simulations are shown in Fig. fT81 
where the transmission coefficient is plotted in depen- 
dence of the chemical potential fi for a fixed repulsive 
nonlinearity g\A\ 2 = 4. Firstly, we notice that the fi- 
nal state ip(x,t en d) is indeed a stationary solution as dis- 
cussed in the previous sections. The transmission coeffi- 
cient |T(^)| 2 follows closely one of the steady state results 
which are computed as described in the previous sections. 
The small deviations are caused by numerical errors, in 
fact mainly by the (still not perfectly absorbing) bound- 
ary conditions. Secondly, the final state tp(x, t e nd) is the 
one with the smallest transmission. This behaviour may 
suppress resonant transport, because a resonance state is 
not populated if different stationary states exist with a 
smaller transmission. This can be seen in Fig. 1181 around 
/X « 5 and [i < 2. However it is possible to access states 
with a higher transmission by a suitable preparation pro- 
cedure as it has been shown in || . If there is no bistabil- 
ity in the vicinity of a resonance (e.g. around /i w 7.5 in 
Fig. I18|) . resonant transport is found. 

A more complicated situation arises when a beak-to- 
beak bifurcation scenario occurs (cf. Fig. |SJ. As an ex- 
ample, we consider the solution of the time-dependent 
NLSE with a source term as in 1|128|> for the same square- 
well potential as above (Vo = —12, a = 4) and a fixed 
chemical potential fi = 2. Again the source strength is in- 
creased adiabatically to a level S(t) = S± so that the 




20 40 60 



FIG. 19: Instability due to a beak-to-beak crossing scenario 
in the transmission through a square- well potential with Vo = 
— 12 and. Shown is the squared modulus of the wave function 
|")/>(a;,t)| 2 for fi = 2 in a grayscale plot (lower panel) and the 
source strength S(t) (upper panel) for comparison (repulsive 
nonlinearity) . 

system assumes a stationary transporting state, where 
g\A\ 2 = 3.5 is well below the critical value for the beak- 
to-beak bifurcation. Then the source strength is again 
slowly increased up to S(t) = S2, so that g\A\ 2 becomes 
larger than the critical value of the beak-to-beak bifur- 
cation. The lowly-transmitting stationary state ceases to 
exist and the time evolution cannot follow the station- 
ary solutions adiabatically any longer. This is shown in 
Fig. El Indeed one observes that the system becomes un- 
stable and finally does not assume a stationary state any 
more. However, the transmission through the square-well 
increases. The beak-to-beak crossing scenario and its im- 
plications on the dynamics will be discussed in detail in 
a subsequent paper psj . 

VIII. CONCLUSION 

The scattering of a BEC by a finite square well poten- 
tial has been discussed in terms of stationary states of the 
nonlinear Schrodinger equation. Neglecting the mean- 
field interaction outside the potential ingoing and outgo- 
ing waves and thus amplitudes of reflection and trans- 
mission can be defined. The transmitted flux shows a 
bistable behaviour in the vicinity of the resonances. For 
repulsive nonlincaritics, the transmission coefficient also 
shows beak-to-beak bifurcations. Wherever these beak- 
to-beak scenarios occur, the scattering phase shows loops. 

An analytical expression for the position of the res- 
onances is derived. The transitions from resonances to 
bound states and vice versa due to attractive/repulsive 
mean-field interaction are proven analytically and an ex- 
plicit formula for the transition interaction strength qt is 
found. Furthermore, the positions of the bound states of 
the system arc calculated and compared with the case of 
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vanishing interaction g = 0. It is found that the oscilla- 
tion theorem also holds for interaction parameters g ^= 0. 

Finally, we solve the time-dependent NLSE numeri- 
cally in order to discuss the dynamics of the system for 
g > 0. If the amplitude \A\ of the incoming flux is slowly 
ramped up for a fixed value of the chemical potential /x, 
the system populates the lowest branch of the stationary 
transmission coefficient, as it was first shown in for 
a double barrier potential, so that resonant transport is 
usually be suppressed. In the vicinity of a beak-to beak 
bifurcation, the system can become unstable if g\A\ 2 be- 
comes larger than the critical value of the beak-to-beak 
bifurcation so that the lowest branch of the stationary 



transmission coefficient suddenly ceases to exist. These 
beak-to-beak crossing scenarios will be further investi- 
gated in a subsequent paper [2jj . 
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